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Abstract 

In a recent paper titled Coherent electromagnetic wavelets and their 
twisting null congruences, I defined the local inertia density I {x,t), 
reactive energy density Ti{x, t), and energy How velocity v{x, t) of an 
electromagnetic field. These are the field equivalents of the mass, rest 
energy, and velocity of a relativistie particle. Thus TZ — Ic^ is Lorentz- 
invariant and \v\ < c, with equality if and only if 7?. = 0. The excep- 
tional fields with \v\ — c were called coherent because their energy 
moves in complete harmony with the field, leaving no inertia or re- 
active energy behind. Generic electromagnetic fields become coherent 
only in the far zone. Elsewhere, their energy fiows at speeds v{x, t) < c, 
a statement that is surprising even to some experts. The purpose of 
this paper is to confirm and clarify this statement by studying the 
local energy flow in several common systems: a time-harmonic elec- 
tric dipole field, a time-dependent electric dipole field, and a standing 
plane wave. For these fields, the energy current (Poynting vector) is 
too weak to carry all of the energy, thus leaving reactive energy in its 
wake. For the time-dependent dipole field, we find that the energy 
can flow both transversally and inwards, back to the source. Neither 
of these phenomena show up in the usual computation of the energy 
transport velocity which considers only averages over one period in the 
time-harmonic case. 
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1 Introduction 

Given an electromagnetic field {E,B) in vacuum, its energy density U and 
energy current (Poynting's vector) S are defined in Heaviside-Lorentz units 
(eo = ^0 = 1) by 

U{x) = ^{E{xf + B{xf) and S{x) = E{x) x B{x), (1) 

where x = {x,t) £ are the spacetime coordinates. By Maxwell's equa- 
tions, the four- vector field {S,U) satisfies Poynting's conservation law 

dtU + cV ■S = -E-J, (2) 

where c is the speed of light and J is the current density. In the absence of 
a current, energy is conserved. Due to the vector identity 

{E X Bf = E^B^ - {E ■ Bf, 

we have 

U^-S^ = ^{E^-B^f + {E-Bf>0. (3) 

At every event x, {S,U) is either future timelike {U{x) > \S{x)\) or future 
lightlike {U{x) = \S{x)\). 



2 



The inequahty (3) does not depend on Maxwell's equations. It is satisfied 
by any two vectors {E,B). To understand its significance, note that the 
electromagnetic momentum density is given by^ S{x)/c, so {S/c,U) is the 
energy-momentum density of the field. For a relativistic particle with energy 
E and momentum p, we have the analog of (3), 

E"^ - c^p^ > 0, 

The mass m and velocity v of the particle are given by 

E/c^ ~ ~E ' 



m = c~^\/ E'^ — c^p^ and v = = (4) 



which implies that v"^ < c^. Thus it makes sense to define [Kll] the electro- 
magnetic inertia density I{x) by 



I = c-Wu^-S^ = ^^{E^-B'^)^ + A{E-B)^ 



(5) 



and the electromagnetic energy Bow velocity v{x) by 



cS 


^ V = \v\ < c. 


V = — = 


U 





(6) 



It will be useful also to define the reactive (rest ) energy density 

TZ{x, t) = c^I{x, t) = ^yU{x, tf - S{x, tf. (7) 

Since E^ — B^ and E B are the two Lorentz invariants of the field, X{x) and 
lZ{x) are local, Lorentz-invariant (scalar) spacetime fields. Their physical 
significance, as well as that of v{x), is the subject of this paper. Our basic 
theme can be summarized by the following local statement in spacetime: 

// |>S'(x)| < U{x), then the energy flow at x does not carry away all of the 
energy, leaving positive rest (reactive) energy and inertia densities at x. 

Reactive energy is an important topic in antenna theory [Y96] . If an antenna 
generates a great deal of reactive energy, this slows down the transmission 
of energy and makes the antenna inefficient. Effectively, IZ forms an elastic 
'cushion' between the source and the far-zone radiation field. While reactive 
energy is defined only for narrow-band fields in the literature, our definition 
(7) extends this concept to arbitrary electromagnetic fields. In fact, since 



^In SI units [J99, page 261], the electromagnetic momentum density is S/(? 
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narrowband approximations are closely related to nonrelativistic limits [K96] 
and TZ is Lorentz-invariant, it can be said to be the ultimate 'wideband' 
definition of reactive energy. 

The quantities S and TZ can be expressed succinctly in terms of the pair 
of complex conjugate vector fields 

F{x) = E{x) + iB{x), F{x)* = E{x)-iB{x) (8) 

as follows: 

U=^-^ S= 7^ = ^ (9) 

2 2i 2 ^ ' 



where 



|_p|2 = F* -F and F^ = F ■ F. 



The combinations E it iB have been called Riemann-Silherstein vectors 
[B3] and Faraday vectors [B99]. They have been rediscovered many times 
and were used extensively by Bateman [B15]. See also [K3], where F{x) is 
continued analytically to complex spacetime, and [K4]. 

By (5), 

n{x) = ^ E{xf - B{xf = E{x) ■ B{x) = = 0. (10) 

An electromagnetic field with F{x)^ = is said to be null at x. Nullity is a 
local, Lorentz-invariant property. It is the field counterpart of masslessness 
in a particle. Indeed, 

v{x) = c ^ X{x) = 0. (11) 

Electromagnetic energy Hows exactly at the speed of light only at events x 
where the field is null. Elsewhere, it Bows at speeds less than c and has a 
positive inertia density. 

Although this simple fact should be widely known in classical electrody- 
namics, I've been unable to find any clear reference to it in the mainstream 
literature and in discussions with several knowledgeable colleagues. The sole 
exception, to my knowledge, is a brief note in [B15, page 6]. 

However, a kind of average energy flow velocity per period of a time-harmonic 
EM fleld is well known in the literature under the name energy transport 
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velocity. ^ It is connected to the above instantaneous flow velocity as follows. 
A time-harmonic field is given by 

E{x,t) = Re {e-''^*E^{x)), B{x,t) = Re {e-''^^B^{x)), (12) 

where E^{x) and B^{x) are complex Fourier components at a positive 
frequency u. Thus 

S{x,t) = S^{x) + S'^{x,t), (13) 

where 

S^{x) = ^{E^ X B*^ + El X B^) = ^ Re {E^ x B*J (14) 
is the average of S{x,t) over one period 2tt/ljj and 

S'Jx, t) = ^Re {e-^'^'E^{x) x B^{x)) (15) 
oscillates at the optical frequency 2lo. Similarly, 

U{x,t)=U^{x)+Ui{x,t) (16) 

where 

U^{x) = -^[\E^{x)\^ + \B^{x)\'') (17) 
is the average of lA over one period and 

Ui{x,t) = ^Re [E^{xf + B^{xf]) (18) 

oscillates at 2a;. 

The energy transport velocity is now defined as the ratio of the averages: 

By comparison, our instantaneous energy flow velocity (6) is given in terms 
of the real fields E, B by 

''^^''^"'zY^(a;)+ZY^(:r,t) -^''E{x,ty + B{x,t)^- 



thank Professor Andrea Alu for pointing this out. 
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Obviously (20) cannot be recovered from (19), but since vibrations at optical 
frequencies are generally unobservable, it might be argued that Vi^{x) suffices 
for all practical purposes. 

However, it cannot be claimed that Vi_j(x) is the time average of v{x,t) 
because the ratio of averages is generally not equal to the average of ratios. 
While is more easily computed than the time average v, it is not a 
good approximation to the latter under all circumstances. 

Furthermore, v^^ ignores some time-dependent aspects of energy transport 
which help explain how the energy can flow at speeds v < c even though the 
waves communicating it propagate at c. This will be illustrated by studying 
several well-known systems. 

2 Time-harmonic electric dipole field 



The field of an oscillating electric dipole [J99] with frequency a; > is given 
by (12) with 

EUx) = k^—e{x), e(a;) = 2A(32-i3)f + (l + i3-32)p (21) 
47rr 

B^{x) = k^—b{x), b{x) = {1 + ii)f X p, 

where k = uj/c is the wavenumber, p is a real electric dipole moment, 
Pj_ = rx{pxr)=p—Xr, \ = f ■ p = p cos 6, p =\p\ 



is a dimensionless zone parameter: 3 — t- oo in the near zone, and 3 — t- in 
the far zone. A straightforward computation gives 



and 
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1 

kr 



e X b* = (1 + if)p^sm^er + 2ii{l + f)p cos 6 p ^ 



(22) 



|e|2 + |fo|2 = p2[2 sin2 e + 432 cos^ 9 + 3^(3 cos^ 9 + 1)]. 
Since all but the first term in e x 5* is imaginary, we get 




(23) 
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The average energy flow velocity over one period is therefore 



cr 



V„AX 



1 + 232 cot2 e + ^3^(3 cot2 e + csc2 I 



r. 



(24) 



The dipole energy thus flows with an average speed v^j < c, approaching c 
only in the far zone. In the near zone it clusters near the equatorial plane 
9 = Tr/2. On any sphere of constant r, vanishes at the poles and increases 
monotonically towards the equator, as shown in Figure 1. 




Figure 1: The average energy flow velocity v^^ per period in the near zone (left, 
3 = 2) and the far zone (right, 3 = .05). v^^ is small and flat in the near zone, 
approaching c in the far zone except for the poles, where it vanishes. 



Equations (23) show that while r'^S^j is zone-independent, r'^Uuj increases 




Figure 2: The average reactive energy density TZ^i with p — k — 1 in the near zone 
(left, 3 = 5) and the far zone (right, 3 = .5). hi the near zone, it clusters around 
the z-axis because w^j = there, so the most energy is left behind. 
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monotonically as we approach the origin. We interpret 

n^{x) = c^X^ix) = ^U^ixY - S^ixf (25) 
as the average reactive energy density per period. 

While the radiating wave front propagates at c, energy is left behind wher- 
ever \Si^{x)\ < h(u){x). The abandoned energy is reactive. 

Plots of TZi^ in the near and far zones are given in Figure 2. 

3 Time-dependent electric dipole field 

This example explains the seeming contradiction that while the time-domain 
dipole fields propagate at speed c, their instantaneous energy density gen- 
erally flows at \v\ < c. Consider a general time-dependent electric dipole 
moment f{t)p fixed at the origin. Its field is 

1 r°° 

E{x,t) = -Re du;e-"^^f{uj)E^{x) 
71" Jo 

1 

Bix,t) = -Re due''"^^f{uj)B^{x) 
^ Jo 

where / is the Fourier transform of /. Inserting (21), we find 

E{x,t) = i{r,tr)pcos6 f — e{r,tr)p± B{x,t) = b{r,tr)p x f (26) 
where tr = t — r/c is the retarded time and 

' ATTC^r A'Kcr'^ 

e(r,t,) = 6(r,t.) + {i^. 

(28) 

Note that i represents the longitudinal electric component of the field while 
e and h represent its transversal electric and magnetic components. Thus 

F = E + iB = £p cos 9r — epj_ + ibpj_ x f (29) 



8 



and we obtain 



S = ebp^ siv? 6f + b£p cos 6 pj_ = SrV + S±pj 



2' ' 2 
S has some remarkable properties 



[{e^ + }?)s\n^e + l^ cos^ 9]. 



(30) 
(31) 



• The second term in (30), determined by the radial component is 
transversal. It corresponds to two of the imaginary terms in (22) 
which were left out when we computed the average energy flow per 
period in (23). Hence not all of the energy is flowing outwards. 

• In the far zone we have 

so the radial term in S dominates and the energy is flowing outwards. 

• In the near zone, eb may be negative, hence energy may flow inwards 
as well as transversally. This will be confirmed below. 

The expression for the instantaneous speed v = c\S\/U is rather complicated. 
However, the reactive energy density 



7^ 



^\{e^-b^)sm^9 + e^cos^^ 



gives a simple expression for v in the relativistic form 



n 

U 



6^) sin^ I 



+ ^2 cos^ ( 



(e2 + 62) sin2 9 + 1"^ cos^ 9 ' 



(32) 



(33) 



Recall that for the time-harmonic dipole we had = Q on the dipole axis. 
This extends to the instantaneous velocity for the pulsed dipole field: 



sin 61 = 



1 



1 



v = 



provided i ^ 0, which can occur only at special values of r and t as shown 
below. Since the expression on the right of (33) is nonnegative, we have 
f < c as required. Furthermore, 



v = c ^ (e^ -b'^)sin^9 + fcos^9 = 0. 
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For a generic pulse /, this can only be satisfied on sets of zero measure in 
spacetime. For example, 

v = c if - 6^ = £ = 0, (34) 

which looks a little like a scalar version of the nullity condition (10). Such 
conditions can be satisfied only for special values of r, and t. For example, 
if £(r, t) = for a Unite time interval, then 

f{t) = --f{t) f{t) = Ae-'^"\ 
r 

This is unacceptable as a pulse function because it depends on r.'^ This 
shows that I can have only isolated zeros. Similarly, — 6^ can have only 
isolated zeros. Thus we have 

v{x,t) < c almost everywhere (a.e.) in spacetime. (35) 

The energy of the pulsed dipole field flows at speeds less than c a.e. 

This statement feels uncomfortable because (26) and (27) show that the 
fields E,B, hence also U and TZ, depend on t only through the retarded 
time t — r/c. If / is a sharp pulse, then so are i, e, b. It follows that U and 
S, like the fields, are small unless r « ct. How can this be reconciled with 
(35)? To investigate this, consider the modulated Gaussian pulse 

f{t) = g{t)cosujt where g{t) = e"''*'/^ k>0, u > (36) 

and set c = 1 and p = 1 for notational convenience. Then 

f'{t) = —Ktg{t) cosut — ojg(t) sinwt 

f"{t) = {k^I^ — k — ijj^)g{t) cosujt + 2KUjtg{t) sinwt, 

hence 

A7:rH{r,t) = ^i,{r,t)g{t) (37) 
W6(r,t) = nb{r,t)g{t) 
iTrr^e{r,t) = ^ieir,t)g{t) 

with 

fii{r,t) = 2(1 — Krt) cos ujt — 2u}rsmujt (38) 
/Uft(r, t) = [K^r'^t^ — (k + a;^)r^ — Krt] cosut + [2KUJr'^t — ur] sinut 
fJ-e{r, t) = lJ,b{r, t) + cos wt. 

^It also grows exponentially for t < 0, but we could set f — for t < 0. 
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Six 


t) = 




So{x 


tr) 


U{x 


t) = 


'giuy 


2 


tr) 


n{x 


t) = 


'gitrY 


2 


,tr) 



(39) 



where 



S q{x^ ^r) 

IAq{x, tf ) 

'^oip^i tr) 



'f + flhUi COS 9pj^ 



1 

2 
1 

2 

with all expressions on the right evaluated at the retarded time 



(40) 



sm 



+ /if cos^ I 



While 



S^IA and TZ all contain the Gaussian factor g{tr) 
in the ratios 

/XeAtb sin^ 9f + fifj/J^e cos^ 6 



this factor cancels 



V 



2c- 



(/^e + sin^ 6* + /i| cos2 i 



and 



1 



^^^) 



sm 



+ cos^ ( 



(a^c + A^b) sin^ 6* + /if cos2 6 



(41) 



(42) 



These ratios depends only on the relative sizes of the factors fJ.e,fJ-b and fig 
at time tr- The pulse g{tr)^ cancels in (42). 

It can be argued that for large It^l, both S and U are extremely small and 
hence their ratio has little meaning. When almost no energy is flowing, the 
speed of its flow is largely of academic interest. This is certainly true in the 
general case (33) at times when TZ and U vanish identically. 

Equation (42) confirms that v vanishes on the dipole axis (sin 6 = 0). It also 
vanishes for special values of r, 6 and t, for example when fif, = or when 
= = 0. There, the instantaneous energy is entirely reactive. 

Remark 1. In the discussion below Equation (30) we have noted that when 
eb < 0, the field energy flows inwards. For the modulated Gaussian pulse 
(36), this occurs when 

Heir, tr)flbir, tr) = fJ-bir, tr)"^ + /^^(r, tr) COSUjtr < 0. 
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Since iJ,b{r,tr) = 0{r'^) as r — oo at given tr, it follows that 

r — )■ OO =^ fJ-eif, tr)fJ-b{r, tr) ~ fJ-bi^, tr)^ > =^ > 0. 

But in the near zone, the radial velocity Vr in (41) can be negative, in 
which case the energy flows inwards. Figure 3 shows the behavior of Vr as 
a function of r on the equatorial plane 9 = tx 12. Note that the energy flows 
inwards for small values of r and then quickly transitions to flowing outward 
at Vr ~ c. The ingoing flow does not show up in the time-averaged velocity 
(24), which is outgoing in all zones. 




Figure 3: Plots of the radial energy flow velocity Vr for the Gaussian pulse (36) as 
a function of r with = (left) and = 1 (right). We have set 9 = tt/2, ui = 1 
and K = 0.1. When v,. < 0, the energy flows towards the origin. 

Remark 2. For a general pulse function f{t), no simple relations like (37) 
exist. Nevertheless, the velocity v = cS /lA depends on the relative sizes 
of f",f' and / through the ratio (33), and not merely on f{t). While the 
pulse does not actually cancel, as it did for the Gaussian, the numerator 
and denominator of (33) are both small when tr is sufficiently large, hence 
nothing definite can be said about v without knowing the relative sizes of 
I, e, b. This explains how the energy can flow at speeds less than c while the 
fields propagate at c. 

However, the value ofv{x,t) says nothing about the quantity of energy flow- 
ing at this velocity. For example, we find that 

sin 6* / and fct^ > 1 =^ v ^ c. 

Although almost no energy remains when tr is large, what little there is 
fiows nearly at the speed of light provided we are not on the dipole axis. 
As an extreme case of this, let /(t) vanish outside the interval T = [ti,t2]. 
Then v{x,t) is undefined for t ^T. 



12 



1 



Figure 4: Left to right: The reactive energy density TZ of the modulated Gaussian 
dipole field with r = 1, w = 1, /t = 0.1 at t = 0.5, 1, 1.3, 1.5, 3. The vertical lobes 
and horizontal tubes are due to the sin^ 6 and cos^ 9 terms in (40), which take turns 
dominating because fii, fj-b and /ig all oscillate with period tt. 



Figure 5: Left to right: The energy flow speed v in (42) of the modulated Gaussian 
dipole field with r — 1, u — 1, k — 0.1 at t = 0.5, 1, 3, 4, 30. Note that v has a zero 
near i = 3, and u — c as i — ^ oo even though TZ,U — 0. 



Figures 5 and 4 show the reactive energy density energy TZ and the energy 
flow speed 



of the modulated Gaussian dipole field on the unit sphere r = 1 at various 
times. As t — )• oo, u — )■ c outside the dipole axis even though the quantity 
of energy flowing vanishes. 

4 Standing plane wave 

Our final example demonstrates an important feature of null fields: they do 
not interfere with themselves as they propagate. 





(43) 
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The simplest examples of null fields are plane waves with real wave vectors.^ 
Consider a pair of linearly polarized plane waves propagating along the 
positive and negative z-axis, 

E± (x) = xE cos(ct =F z) , B± = ±yE cos(ct =F z) . 

These are solutions of Maxwell's equations with complex representations 

F± = E± + iB± = {x±iy)Ecos{ct^ z). (44) 

The energy-momentum of F± is 

U± = E"^ cos^ {ct^z) S± = ±zE'^ cos^ (ct =F 2) , 

hence the inertia density and energy flow velocity are 



^ 1 , „9 , cS-i- 



Thus F± are null fields moving in the zizz direction at speed c. 
Now consider the standing wave 

F = F^ + F^ = {x + iy)Ecos{ct - z) + {x - iy)Ecos{ct + z), (45) 

which can also be written as 

F = E + iB where E = 2xE cos ct cos z, B = 2yE sin ct sin z. 

Note that (45) is expressed in terms of the traveling waves cos(ct it z) while 
E and B are not. It will therefore be more revealing to use F. We find 

I(z,t) = ^ = 2(E'^/c'^)\cos(ct- z)cos(ct + z)\ 
F* ■ F 

U{z, t) = — - — = E^ [cos^ {ct - z) + cos^(ct + z)] 



2 

F* X F 

~2A 



S(z, t) = — — — = zE"^ rcos^(ct - z) - cos^(ct + z) 



*In the plane-wave spectrum representation [HY99] (also called the angular spectrum 
representation) , electromagnetic fields are expressed as superpositions of plane waves with 
complex wave vectors. Such 'inhomogeneous' plane waves are not null since they include 
non-propagating evanescent waves. 
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Figure 6: Plots of v{z,t) at ti = n/ld {left) and half a cycle later, at t2 = 97r/16 
{right). On the left, the energy travels to the right in the intervals (0,7r/2) and 
(7r,37r/4), and to the left in {tt/2,tt) and (37r/4, 27r). On the right, it travels in 
the opposite directions. Furthermore, v{z,0) — v{z,tt) = and v{z,t) jumps 
discontinuously across the odd nodes Z2n+i as i — and across the even nodes Z2„ 
as t — > TT. Note that v is not at even close to harmonic. 



and the energy flow velocity is 



v{z,t) 



cS 

U 



zv{z,t), v{z,t) 



cos^(ci — z) — cos^(ct + z) 
cos^(ci — z) + cos^(ct + z) 



(46) 



As expected, \v\ < c. Note that v is periodic of period vr/c, and 

V = <^ cos^ {ct — z) = cos^ {ct + z) <^ ct + z = ±{ct — z) + mr. 

Hence v has fixed nodes in space and time, as seen in Figure 6: 

nvr mr 
v{z,t) = 4^ z = — = Zn or ct = — = ctn, n G Z. (47) 

Since v{z, t) changes sign at z„, the energy is reflected at the nodes. 

The energy oscillates back and forth between successive nodal planes, and 
v{z,t) oscillates between ±c at any z which is not a node. 

For non-nodal z, we have 

cos^(ct + z) = 4^ z {t) = {n + l)7r - ct 4^ v = c (48) 
cos^(ct - z) = ^ z{t) = {n + ^)'K + ct ^ V = -c. 

That is, if we begin at the node z(0) = Z2n+i = (n + ^)7r and travel to 
the left at speed c, we see the velocity of -F+ except when crossing the 
nodes. If instead we travel to the right at speed c, we see the velocity of 
F-. This behavior can be seen in Figures 6 and 7. Note that v{z,t) jumps 
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4 2 



Figure 7: Plot of v{z,t) (46) showing the nodes (47) and confirming (48). 



discontinuously between itc across the nodes. This is because rather than 
crossing the nodes, the velocity builds up to itc and is then immediately 
reflected to =Fc. 

Had we used the time-averaged formula, we would have obtained the ex- 
pected result = 0. However, (6) gives a detailed, instantaneous picture 
of the movement of energy. 



5 The incoherence of electromagnetic fields 

The above examples show that we must distinguish between the propagation 
speed of the fields {E, B) and the flow velocity of their energy. Let us think 
about this from a purely mathematical point of view. 

• The fields satisfy homogeneous wave equations outside their sources, hence 
they propagates at speed c. 

• On the other hand, {S,U) are quadratic functions of the fields and satisfy 
Poynting's conservation law outside of sources, 

dtU + cV ■ S = dtU + V ■ (vU) = 0. (49) 

Hence U behaves like the density of a compressible Euid flowing with velocity 
V. Although (49) follows from Maxwell's equations, there is no a priori 
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reason why the two speeds should be equal. Indeed, as we have shown, they 
are in general different; only in the far zone^ does w — t- c. 

A rough way to understand why f < c is by analogy with water waves. The 
mass carried by the waves has a definite speed at each point and time, but 
this need not coincide the propagation speed of the wavefronts. 

But null fields are an exception: all of the energy is carried away by the wave 
since v = c. For them, the energy flows in unison with the waves, hence they 
play a very special role in electrodynamics. 

The following definition is motivated by the fact that the electromagnetic 
energy propagates coJierentiy with the field if and only if F{x)^ = 0. In fact, 
the traveling plane waves (44) are mutually coherent in this sense even at 
different spacetime points: 

F±{x) = {x±iy)Ecos{ct^ z) F±(x) • = Vx, y G 

since {x it iy)"^ = 0. 

Definition 1 The incoherence of an electromagnetic field F = E + iB is 
the complex function 



3{x, y) = E{x) ■ E{y) - B{x) ■ B{y) + iE{x) ■ B{y) + iB{x) ■ E{y). 



3 measures the incoherence of F across space and time by comparing F at 
two different events. It can be related to the coherence functions of statistical 
optics [W7] as follows. 

If F is time-harmonic. 



then the average 3u] of the equal-time incoherence function over one period 
is given by 



y) = e{x) ■ e{yy + e{xY ■ e{y) - b{x) • &(y)* - b{xy • b{y) (51) 
+ ie{x) ■ b{y)* + ie{xY ■ b{y) + ib{x) ■ e{y)* + ib{x)* ■ e{y). 



^Provided, of course, that a far zone exists. When the sources have an infinite extent 
or are 'at infinity' as in the case of plane waves, no 'far zone' may exist. For example, 
the standing wave (45) has no far zone, and the traveling plane waves (44) have all of 
spacetime as their far zone. 




(50) 



which is expressed in terms of the fields E,B as 



E{x,t) = 2Re (e~"^*e(^)), B{x,t) = 2Re (e"*'^*6(iE)) 
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This can be expressed in the compact form 



(52) 



where 



f _^_{x) = e{x) + ib{x) and f _{x) = e{x) — ib{x) 



are independent because e and b are complex. 

The connection with the coherence functions of statistical optics is obtained 
by extending J to random fields, where we have an ensemble of fields F(x) 
and its incoherence is defined by the ensemble average 



Equation (52) shows that for a random time-harmonic field, ( 3 is a specific 
combination of coherence functions for the electric and magnetic fields. But 
while coherence functions are designed to test the correlation of a set of fields, 
our incoherence function represents their electric-magnetic imbalance, in the 
sense that a random field F is perfectly balanced at x when 



{F{xf) = 0, i.e., {E{xf) = {B{xf) and {E{x) ■ B{x)) = 



This notion of balance is required for the field's energy to propagate co- 
herently with the field. In the time-harmonic case, this means we have the 
following identities between the correlation functions: 



Whereas the coherence functions of statistical optics are defined only for 
random fields, we have seen that the incoherence function has a deep signif- 
icance even for a deterministic field. 

A generic electromagnetic field in free space is null along a set of 2-dimensional 
hypersurfaces S in spacetime since F{x,tf = imposes two real conditions 
on the four spacetime variables {x,t).^ The time slices St of S are (generi- 
cally) curves in space, and these curves evolve with t.'' Thus, when a field 

®We have seen an example of this in Section 3, where the coherence condition (34) for 
the time-dependent electric dipole field can only be satisfied on isolated surfaces. 

'^Bialynicki-Birula [B3] calls these moving curves electromagnetic vortices, but the 'ro- 
tations' around these curves are duality rotations F — >■ e"''f rather than rotations in 
physical space. 




(53) 



Re{e{x)-e{yy)= Re {b{x) ■ b{yy ) 
Re ( e{x) ■ b{yy ) = - Re ( b{x) ■ e{yy ). 
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is not null in an extended region of spacetime, its energy flows at the speed 
of light only along such curves. Elsewhere it flows at speeds v < c, although 
f — )• c in the far zone. 

So far, our only example of a null field has been the traveling plane waves 
F± (44). It is easy to make a plane wave null because its energy flow 
velocity is a constant vector. Various other globally null electromagnetic 
fields are known where F'^ vanishes almost everywhere, with the possible 
exception of singularities whose supports have zero measure. Such fields 
play an important role in the solution of the Einstein-Maxwell equations of 
general relativity [B15, R61, T62, RT64]. However, no extended, compactly 
supported sources appear to be known which radiate null fields everywhere 
outside the source region. Such null field, called coherent electromagnetic 
wavelets, was recently constructed [Kll]. It is radiated by a relativistically 
spinning charged disk, and the radiation follows a twisting null congruence 
of light rays. This is a space-filling set of null lines (world lines of particles 
traveling at speed c) which are determined by the energy fiow velocity vector 
V = cS /U with \v\ = c The existence of such null congruences is the key 
to the construction of nontrivial null fields. 
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